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Abstract : Compact Heat Exchangers (CHEs) are 
increasingly being used on small and medium scale 
industries. Due to their compact size and efficient design, 
they facilitate more efficient heat transfer. Better heat 
transfer would imply lesser fuel consumption for the 
operations of the plant, giving improvement to overall 
efficiency. This reduction in consumption of fuel is a step 
towards sustainable development. This report exclusively 
deals with the study the spiral heat exchanger.The design 
considerations for spiral heat exchanger is that the flow 
within the spiral has been assumed as flow through a duct 
and by using Shah London empirical equation for Nusselt 
number design parameters are further optimized.This is 
accompanied by a detailed energy balance to generate a 
concise mathematical model. 

1) Introduction 

Heat Exchanger is an apparatus that has two streams of fluids 
at different temperatures that exchange heat in order to satisfy 
the heating or cooling requirements of the system They are 
installed in plants to help improve energy utilization. The 
energy from the streams at high temperature is transferred to 
streams at low temperature whenever required. Heat transfer 
only occurs when the system satisfies both the first and second 
law of thermodynamics. 

A spiral heat exchanger (SHE), may refer to a helical (coiled) 
tube configuration, more generally, the term refers to a pair of 
flat surfaces that are coiled to form the two channels in a 
counter-flow arrangement. Each of the two channels has one 
long curved path. A pair of fluid ports is connected 
tangentially to the outer arms of the spiral, and axial ports are 
common, but optional. 

The concept of a spiral heat exchanger is as simple as it is 
sophisticated. Two or four long metal strips, onto which spacer 
studs are welded, are wound around a core, thus creating two 
or four equally spaced single passage channels. The concentric 
shape of the flow passages and the studs yield turbulence 
already at low Reynolds numbers. By optimizing the flow 
pattern heat transfer is enhanced, whilst fouling is reduced. 

The SHE is its highly efficient use of space. This attribute is 
often leveraged and partially reallocated to gain other 
improvements in performance, according to well-known 
tradeoffs in heat exchanger design. A compact SHE may be 
used to have a smaller footprint and thus lower all-around 
capital costs, or an over-sized SHE may be used to have less 
pressure drop, less pumping energy, higher thermal efficiency, 
and lower energy costs. 



The most fascinating feature of SHE is that due to its 
geometry, there is very limited fouling. As the fluid flows in a 
spiral, deposition on the walls f the heat exchanger is 
negligible. This means lower down time as compared to the 
conventional Shell and Tube heat exchangers. 




(Fig 2.1) conceptual flow of fluids (hot and cold) in SHE, 
courtesy Alfa Laval) 

While designing any heat exchanger for a plant, the heat flux, 
temperature drop and the surface area required to transfer the 
heat is known via the preliminary Heat Exchanger Network 
calculations. The design parameters aim to achieve these 
requirements in the most optimal way plausible keeping in 
mind safety, operability, maintainability, sustainability and 
profitability. In this paper, I assume that the heat flux and the 
inlet-outlet temperature of both hot and cold streams are 
known. 

2) Methodology 

For the design equation calculation, a differential cross section 
is considered for the spiral heat exchanger and then used 
counter current flow LMTD method to solve for the design 
parameters. 

/th*t^t--*\ (Th.out — T c i n )— (Th j n — T c out ) 

( LMTD ) = , ^h.out-Tc.in) , 

( r h,in _T c,out J 

T hin = Hot Stream Inlet Temperature 

T h out = Hot Stream Outlet Temperature 

T c in = Cold Stream Inlet Temperature 

T c ,out - Cold Stream Outlet Temperature 

q = UA {LMTD) (Eql.2) 
Here, the heat flux is represented as a function LMTD. 
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Mathematical Modeling Spiral Heat Exchanger 
Design: The distance between the sheets in the spiral channels 
are maintained by using spacer studs that were welded prior to 
rolling. Once the main spiral pack has been rolled, alternate 
top and bottom edges are welded and each end closed by a 
gasketed flat or conical cover bolted to the body. This ensures 
no mixing of the two fluids occurs. 



Assuming the flow is through a duct at each location. The 
resistances (as represented in fig. 2.3) to the heat flux are 
offered by the thermal resistivity of the separating wall and 
convective heat flux offered by the flowing fluids. As the heat 
flux for every compartment flows from the hot fluid to the 
cold fluid, the overall heat transfer coefficient can be evaluated 
using the concept of resistances in series. Estimating the over 
all heat transfer coefficient using this concept. 







HOT FLUID 



COLD FLUID 



2a 



2b (Fig 
2.2) CROSS SECTION CONSIDERING HOT, COLD 
ELEMENTS EACH. 
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(Fig 2.3) Block diagramof thermal resistances in series 

(Eq2.6) 



X K. h/ 



The overall heat transfer coefficient 



Using LMTD (log mean temperature difference) Method [4] 

Assuming counter current flow of fluids 

Heat balance equation (taking LMTD approach) 

dq = —rri h C h dT h = m c C c dT c = U(j h — T c )dA 

(Eq2.1) 
l l 

dTh -dT c = -dq{— +— ) 
d(T h -T c )=-U(.T h -T c X-^+-^r)dA 

J1 (T h -T c ) J ° \m c C c m h C h ) 



ln[% r) - -UA ( 1- ) 

yjh,in~ T c,out) m c C c m h C h 

As, 



rri h C h = 



We get, 



T h,in T h,<mt 



' T cin 



1 1 . (T h ,out T c i n 

)-(Th,in-T c ,out) 

q = UA (TT — : -t . -i 

1 , tv A' h.out 1 c,m> . 

hi! fe z 



(Eq2.2) 

(Eq2.3) 
(Eq2.4) 



In Eq2.6, w is fixed as a design parameter and has fixed values 
based on material of construction. K s v& fixed as well as it 
depends on the material of construction which is a design 
decision. ft e andft c both are dependent on the flow of the fluid, 
Nusselt number ( a dimensionless number) is used for the 
estimation of ft (convective heat transfer coefficient). The 
focus here is to estimate overall heat transfer coefficient by 
relating it to the Nusselt number. 



Nusselt number, Nu 



h.l e f f ective 



Kftuid 

(Definition of Nusselt number) 

* fluid 



=>h = 

^effective 

Now, estimating each term on the R.H.S. 



(Eq2.7) 



K fluid i i s a property of the fluid flowing in each cold and hot 
side, it has no control on the design parameters. 

Ineffective estimation; as the flow of the fluid at each location 
can be assumed to be a flow through a rectangular duct 

_ 44 _ 4(2a*2fe) _ 
L effective ~ ~ ~ 2 ( 2a+2b) ~ D h 



As, a «b 

^effective = * a 



(Eq2.8) 



Where, 
(LMTD) : 



q = UA (LMTD) 



(Th ,out ^c^n )-(Th,in-T c ,out) 

l r r ^h.OUt ~T C.in j ^ 

ln *7T^ „ 



[3] 



(Eq2.5) 



(* h,in~ 1 c.out ) 

Where U can be estimated by considering the thermal 
resistances offered by the heat exchanger. 

Assuming at each location around the spiral, the hot fluid 
dissipates equal heat to the outer and inner (radial 
direction). Effectively, each location dissipates energy to a 
location radially outward from it. 



Using Shah-London empirical equation 
Nu Dh = 8. 235(1 - 2.0421c + 3.0853 c 2 - 2.4765 £ 3 + 1.0578 



c 4 - 0.1861 £ 5 ) 

Where e = - 

b 

Using (Eq2.8)in(Eq2.7) 

h = (Nu Dk * K fluid )/4a 
For hot fluid 

K = ( Nu D k * K hot, fluid )/ 4a 

For cold fluid 

K = ( Nu D k * K cold, fluid )/ 4a 



(Eq2.9) 



(Eq2.10) 
(Eq2.ll) 
(Eq2.12) 
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Putting (Eq2.11)and (Eq2.12) in (Eq2.6) 

U = ( + ^ + 

(Wu Dft * K hotfluld ) K s (Nu Dh *K cold fluid ) 

Putting in (Eq2.5) 



q = A- 



LMTD 



(Eq2.13) 



(Nu Dh * K hotf[uid ) K s (Nu Dh *K coldfluid ) 



Nu Dh , is given by the Shah -London empirical equation, and is 
same for both hot and cold fluid as the size of the duct is same 
for both. 

Here K h ,j U id , Kc,suid , K s , Q ,Ch -> are properties and are fixed 
as a design point of view they are dependent only on the bulk 
temperature of the fluid. Assuming these parameters do not 
change much within the working temperature range on the 
heat exchanger they can be assumed to be constant. 
A, is the surface area of the heat exchanger. 

Higher the value of A, more will be the log mean temperature 
difference and hence better heat transfer but larger A means 
more construction material, and hence a higher capital cost 



Optimizing e parameter 
q a Nu Dh 



(Eq2.14) 



As, Nu D = /(e) = 8.235(1 - 2.0421e + 3.0853 <l - 2.4765 



<-D h 

e 3 + 1.0578 e 4 - 0.1861 e 5 ) 



Maximizing by ignoring higher order terms 

SNu D 



(Eq2.15) 



de 



■= (-2.0421+6. 1706e - 7.4295 e )*8.235 



(Eq2.16) 



Forextrema 

— (-2.0421+6. 1706e - 7.4295 e 2 )*8.235 = 

de 

Hence 
Ase > 

Solving equation (Eq2. 15) 
e « 0.4153 



de 



= 8.235(6.1706 - 14.859e) 



At e= 0.4153 

dNu D 
d ^ 

— — = -2.822 * 10 -3 

a e 

As 

dNu n 
3— A 

— ^— = -2.822 * 10" 3 <0 



So, Nu Dh is maximum ate = 0.4153 

From Eq. 2.15 

Nu Dh (0.4153) = 4.172 

Ande = -=0.4153 

b 

=>a = 0.4153ft 



(Eq2.17) 



(Eq2.18) 



These design parameters should be used to get the optimum 
heat transfer.The optimumdesign equation is : 
Substituting Eq2.18 in Eq2. 13 



<7 = A 



LMTD 



0.3983i) 



K hot,fluid K co Id, fluid) K s 



(Eq2.19) 



Now, estimating the A for spiral heat exchangers. 
A = 2(2b)l=4bl (Eq2.20) 
Estimating length of the spiral: 

Using logarithmic spiral equation for estimating the length 
of the spiral [17] 

For a spiral the locus in polar coordinates is: 

r =oc e pe (Eq2.21) 

for the spiral: 

r = 2a ; at# = 2n 

r = 4a ; at# = An 

Hence, solving for oc and /? we get, 

r=ae2n" (Eq2.22) 
Taking a line integral to estimate the value for 1 



r 2 + A 2 d6 



o 

Where, n is decided by the designer. 

On substituting Eq2.22 and solving the integral 



I 



2an 
ta(2) 



(Eq2.23) 



Substituting in Eq2.20; 



A = 4bl = 



8ban 



ln(2) In (2) 

1 + (—y (e n 2 - i) 



As, a = 0.41536 



= >A = !W« | 1 + ( !^ )2 (e -*fi _ 1} 

In (2) aJ v 2tt j v j 

Substituting in Eq2. 19; 



3.3224 b 2 n 
In (2) 



0.3983ft 



( 1 I 1 \+— 
\ K hot,fluid K cold,fluidJ K s 



(Eq2.24) 



In the above equation q and LMTD are fixed as mentioned 
before. The only variable is b which can be estimated from the 
above equation. Hence, the entire heat exchanger can be 
designed. 

4) Results and Conclusion 
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After studying the heat balance equation and applying the 
LMTD method, the optimum design parameters for a spiral 
heat exchanger were achieved, the final design equation which 
is given as 

/3.3224 b 2 n I ln(2V 7 , n ' n(2) „'\ F .. 



hot.fluid K cold,fluidJ K s 



This paper dealt with the design of a compact heat exchanger, 
heat balance was carried out for the spiral heat exchanger 
taking into account both the conductive and convective modes 
of heat transfer. The design equation has been represented with 
physical parameters that decide the dimensions of the entire 
heat exchanger. The initial heat exchanger network 
calculations of any plant provides the heat flux and the change 
in temperatures of the streams, upon feeding those values to 
the above modeled equations, the design parameters for the 
heat exchangers can be estimated. 

This model helps estimate the width of the spiral 'b'. While 
developing the equation using the Shah London's equation for 
Nusselt number which is dependent on 'e' has been optimized 
to produce a design that would provide the best heat transfer 
while using minimum material to construct the SHE. By this 
approach, a more efficient design can be approached. 
Reducing the material would have direct impact on the capital 
expense of the heat exchanger. 

NOMENCLATURE 

rri h ->Mass flow rate of hot fluid (kg/h) 
C ft ->Specific heat capacity of hot fluid (kJ/Kg.K) 
7)j->Temperature of hot fluid (K) 
m c ->Mass flow rate of cold fluid (kg/h) 
C c ->Specific heat capacity of cold fluid (kJ/Kg.K) 

T c ->Temperature of cold fluid (K) 

t/->Overallheat transfer coefficient (W/m 2 k) 

A->Total internal surface area of the heat exchanger (m 2 ) 

hh->Convective heat transfer coefficient for hot fluid (W/m 2 k) 

h c ->Convective heat transfer coefficient for cold fluid (W/m 2 k) 

w->Width of the separating wall (m 2 ) 

K s ->Thermal conductivity of wall (W/mK) 

Kfluid->Thermal conductivity of the fluid (W/mK) 

Nu Dh -> Nusselt number for hydraulic diameter D h 

Kh,suid->Thermal conductivity of hot fluid. (W/mK) 

Kc^uid^Thermal conductivity of cold fluid(W/mK). 

D h -> Hydraulic Diameter (m) 

1 -> length of spiral (m) 

n -> Number turns of the spiral 
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